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By combining the Minkowski inequality and the quantum Chernoff bound, we derive easy-to- 
compute upper bounds for the error probability affecting the optimal discrimination of Gaussian 
states. In particular, these bounds are useful when the Gaussian states are unitarily inequivalent, 
i.e., they differ in their symplectic invariants. 
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I. INTRODUCTION 

One of the central problems in statistical decision the- 
ory is the discrimination between two different proba- 
bility distributions, intended as potential candidates for 
describing the values of a stochastic variable. In general, 
this statistical discrimination is affected by a minimal 
error probability pW , which decreases with the number 
N of (independent) observations of the random variable. 
The general problem of determining pW was faced by 
H. Chernoff in 1952 [lj. Remarkably, he derived an up- 
per bound, today known as "Chernoff bound", having 
the non trivial property of providing p W in the limit of 
infinite observations (i.e., for N — > +00). Very recently, 
a quantum version of this bound has been considered 
in Refs. @, Hj]. Such a "quantum Chernoff bound" al- 
lows estimation of the minimal error probability p( N > 
which affects a corresponding quantum problem, known 
as quantum state discrimination. In this problem, a 
tester aims to distinguish between two possible quantum 
states of a system, supposing that N identical copies of 
the system are available for a generalized quantum mea- 
surement. The problem of quantum state discrimination 
is fundamental in several areas of quantum information 
(e.g., quantum cryptography Q) and, in particular, for 
continuous variable quantum information Contin- 
uous variable (CV) systems are quantum systems with 
infinite-dimensional Hilbert spaces like, for instance, the 
bosonic modes of a radiation field. In particular, bosonic 
modes with Gaussian statistics, i.e., in Gaussian states 
@, are today extremely important, thanks to their exper- 
imental accessibility and the relative simplicity of their 
mathematical description. 

In the CV framework, the quantum discrimination of 
Gaussian states can be seen as a central task. Such a 
problem was first considered in Ref . , where a formula 
for the quantum Chernoff bound has already been de- 
rived. In our paper, we recast this formula by making 
explicit its dependence on the symplectic spectra of the 
involved Gaussian states. The computational difficulty 
of this formula relies on the fact that, besides the sym- 
plectic spectra (easy to compute), one must also calculate 
the symplectic transformations that diagonalize the cor- 
responding correlation matrices. The derivation of these 



symplectic transformations can be in fact very hard, es- 
pecially when many bosonic modes are involved in the 
process. In order to simplify this computational prob- 
lem, here we resort to standard algebraic inequalities, 
i.e., the Minkowski inequality and the Young's inequal- 
ity. Thanks to these inequalities, we can manipulate the 
formula of the quantum Chernoff bound and derive much 
simpler upper bounds for the discrimination of Gaussian 
states. These bounds, that we call Minkowski bound and 
Young bound, are much easier to compute since they de- 
pend on the symplectic spectra only. Notice that, be- 
cause of this simplification, these bounds are inevitably 
weaker than the quantum Chernoff bound. In particular, 
they are useful when the Gaussian states are unitarily in- 
equivalent, i.e., not connected by any unitary transforma- 
tion (e.g., displacement, rotation or squeezing). On the 
one hand, this is surely a restriction for the general ap- 
plication of our results. On the other hand, inequivalent 
Gaussian states arise in many physical situations, and 
easy-to-compute upper bounds can represent the unique 
feasible solution when the number of modes is very high. 

The structure of the paper is the following. In Sec. |TT] 
we review some of the basic notions about Gaussian 
states, with a special regard for their normal mode de- 
composition. In Sec. IIIII we review the quantum Cher- 
noff bound and re-formulate the corresponding expres- 
sion for Gaussian states. The subsequent Sec. |IV] con- 
tains the central results of this paper. Here, we de- 
rive the computable bounds for discriminating Gaussian 
states by combining the quantum Chernoff bound with 
the Minkowski determinant inequality and the Young's 
inequality. We also provide a simple example in order to 
compare the various bounds. Sec. [V] is for conclusions. 



II. GAUSSIAN STATES IN A NUTSHELL 

Let us consider a bosonic system of n modes. This 
quantum system is described by a tensor product Hilbert 
space Ti^ n and a vector of quadrature operators x T := 
(qi,Pi, . . . ,q n ,p n ) satisfying the commutation relations 



[xi,x m ] = 2iSl im (1 < l,m < 2n) , 



(1) 



2 



where 



k=l 



1 
-1 



(2) 



defines a symplectic form. An arbitrary state of the sys- 
tem is characterized by a density operator p S T>(TL® n ) 
or, equivalently, by a Wigner representation. In fact, by 
introducing the Weyl operator [8j 



(3) 



an arbitrary p is equivalent to a Wigner characteristic 
function 



:= Tr 



or to a Wigner function 



W(x) 



(2tt) 2 



exp (-ix T £) x{£) 



(4) 



(5) 



In the previous Eq. (J5J) the continuous variables x T := 
(qi,Pi, . . . , q n ,Pn) are the eigenvalues of x T . They span 
the real symplectic space JC :— (R 2 ",J1) which is called 
the phase space. 

By definition, a bosonic state p is called Gaussian if the 
corresponding Wigner representation (x or W) is Gaus- 
sian, i.e., 



x{£) = cx p 

W(x) 



1 T 

"2« V « 



ix T £ 



expf-Kx-x^V-^x-x)] 
(27r)"VdctV 



(6) 
(7) 



In such a case, the state p is fully characterized by its dis- 
placement x : = Tr(xp) and its correlation matrix (CM) 
V, with entries 



Vim ■= |Tr[{Ai;, Ax m }p] 



(8) 



where Axi := i/ — Tr(i;p) and {, } is the anticommutator. 
The CM is a 2n x 2n, real and symmetric matrix which 
must satisfy the uncertainty principle 



V + in > 



(9) 



directly coming from Eq. (jTJ and implying V > 0. 

Fundamental properties of the Gaussian states can be 
easily expressed via the symplectic manipulation of their 
CM's. By definition, a matrix S is called symplectic when 
it preserves the symplectic form of Eq. [[2"]). i.e., 

sns T = n . (10) 

Then, according to the Williamson's theorem, for every 
CM V there exists a symplectic matrix S such that 



V = S 



s 1 = s 



(J) Vkh 



(ii) 



where the set {vy, • • • , v n } is called symplectic spectrum 
@. In particular, this spectrum satisfies 



J| u k = VdetV 



(12) 



k=l 



since det S = 1 . By applying the symplectic diagonaliza- 
tion of Eq. (fTTjl to Eq. ([9|) , one can write the uncertainty 
principle in the simpler form (Iol | 



Uk > 1 and V > 



(13) 



A. Normal mode decomposition of Gaussian states 
and its application to power states 



An affine symplectic transformation 
(x, S) : x — * Sx + x , 



(14) 



acting on the phase space JC := (R 2n , fi) results in a sim- 
ple congruence transformation V — > SVS at the level 
of the CM. In the space of density operators T){H® n ), 
the transformation of Eq. (fl4|) corresponds instead to the 
transformation 



P U RlS p Ul s 



(15) 



where the unitary U^.s = D(x)Us is determined by the 
affine pair (x, S) and preserves the Gaussian character 
of the state (Gaussian unitary). As a consequence, the 
symplectic diagonalization of Eq. (fTTjl corresponds to a 
normal mode decomposition of the Gaussian state 



p=U : 



x,S 



Lfc=i 



m 



s > 



(16) 



where 



cr(u k ) 



Uk 



3=0 



Uk 



Vk 



1 



\j) k (j\ (17) 



is a thermal state with mean photon number n k — (y k — 
l)/2 ({|j) fc }°^ are the number states for the fcth mode). 
Thanks to the normal mode decomposition (5c,S, {u k }) 
of Eq. (|16p. one can easily compute every positive power 
of an n-mode Gaussian state p. In fact, let us introduce 
the two basic functions 



:= (x + If ± (x - If 



(18) 



which are nonnegative for every x > 1 and p > 0. Let us 
also construct 



and 



G p (x) 



A p {x) 



2" 



2" 



$"(x) (x + If - (x - If 
<W _ (x + lf + (x-lf 



$p(x) (x + lf~(x 
Then, we have the following 



If 



(19) 



(20) 
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Lemma 1 An arbitrary positive power p p of an n-mode 
Gaussian state p can be written as 



where 



and 



p"=(TrpO)p(p) 



Trp» = Y[ G p {v k ) , 
fc=i 



Proof. By setting 

l + Vk 

1 - rjk 



(21) 



(22) 



(23) 



(24) 



into Eq. (1 1 7|) , we have the following equivalent expression 
for the thermal state 



o-fe) = (!- Vk)J2v 3 k \j)k 01 



By iterating Eq. ([16]) we get 

/^ = r x . s ^(X)K%)fl^, 




(25) 



(26) 



for every p > 0. Then, from Eq. (|25[) . we get 



3=0 



[a(r, k )r = (1 - VkTj2Wk) J \j)k 01 = 



(27) 



which, inserted into Eq. ([26]) . leads to the expression 



U, 



Now, from Eq. (|28p we have 



k=l 



and, applying Eq. (124]) , we get 



Trpf = J] 



2'-' 



L x (^ + 1) P - {u k - If 



(28) 



(29) 



(30) 



which is equivalent to Eqs. (|2"2"|) and ([T9"]) . Then, we can 
easily derive the symplectic eigenvalue u kjP of the thermal 
state o-(rj^) which is present in Eq. ([28]) . In fact, by using 
Eq. J2J), we get 

_\ + rg K + l^ + K-1)^ 

" fe < p ~ 1^ - RTIM^ - A ^ fe) ' (31) 



i.e., fk,p is connected to the original eigenvalue v k by the 
A- function of Eq. ([20]) . Finally, by inserting all the pre- 
vious results into Eq. ([28]) we get the formula of Eq. (f2"T]l . 



Notice that, thanks to the formula of Eq. ([21]) . the 
unnormalized power state pP is simply expressed in terms 
of the symplectic spectrum {vk\ and the affine pair (x, S) 
decomposing the original Gaussian state p according to 
Eq. (p~6|) . In particular, the CM V(p) of the normalized 
power state p{p) is simply related to the CM V = V(l) 
of the original state p = p(l) by the relation 



V(p) = S 



0A p (z/ fc )I fc 



Lfc=l 



(32) 



III. QUANTUM CHERNOFF BOUND 

Let us review the general problem of quantum state 
discrimination (which we specialize to Gaussian states 
of bosonic modes from the next subsection IIII Ap . This 
problem consists in distinguishing between two possible 
states pa a nd p b, which are equiprobable for a quan- 
tum system [111 ]. In this discrimination, we suppose that 
N identical copies of the quantum system are available 
for a generalized quantum measurement, i.e., a positive 
operator valued measure (POVM) [l^]. In other words, 
we apply a POVM to N copies of the quantum system 
in order to choose between two equiprobable hypotheses 
about its global state p N , i.e., 



H A : p 



N 



PA 



PA '-- 



N 



and 



H B :p N = p B 



' Pb 



N 

Pa 



P N b 



(33) 



(34) 



A' 



In order to achieve an optimal discrimination, it is suffi- 
cient to consider a dichotomic POVM {Ea,Eb}, whose 
Kraus operators Ea and Eb are associated to the hy- 
potheses Ha and Hb, respectively. By performing such 
a dichotomic POVM {Ea,Eb}, we get a correct answer 
up to an error probability 



Pj£? = \P{H A \H B ) + \PiHs\H, 
= |Tr(^)+lTr(^ 



(35) 



Clearly, the optimal POVM will be the one minimizing 
Now, since Ea = I — Eb, we can introduce the 



Helstrom matrix 11 



1--=P N B 



Pa 



and write 



pW = |-lTrM. 



(36) 



(37) 
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The error probability of Eq. ([37)1 can be now minimized 
over the Kraus operator Eg only. Since Tr (7) = 0, the 
Helstrom matrix 7 has both positive and negative eigen- 
values. As a consequence, the optimal Eg is the projec- 
tor onto the positive part 7+ of 7 (i.e., the projector onto 
the subspace spanned by the eigenstates of 7 with posi- 
tive eigenvalues). By assuming this optimal operator, we 
have 



Tr 



where 



( 7 £ B )=Tr( 7+ ) ^hlli 



I7II1 : = Tr l7l = TrVrh' 



(38) 



(39) 



is the trace norm of the Helstrom matrix 7. Thus, the 
minimal error probability pW := min Perr is equal to 



P(N) = H 1 - 5 \\Pb 



(XL) 



(40) 



The computation of the trace norm in Eq. (|40|) is rather 
difficult. Luckily, one can always resort to the quantum 
Chernoff bound [3J] 



pw < p£» 



where 



and [3] 



3 W 
QC 



log 



o inf iT r(p>-0 



More simply, this bound can be written as 

-1 N 



P, 



(JV) 



QC 



inf Q s 

0<s<l 



where 



Q s := Tr (p^ 



l-s\ 



(41) 



(42) 



(43) 



(44) 



(45) 



Notice that the quantum Chernoff bound involves a mini- 
mization in the variable s. By setting s — 1/2 in Eq. (|4"4")) . 
one can also define the quantum version of the Bhat- 
tacharyya bound [HI 



P B N) ~\[Tr (^pX^pb-)} 



N 



which clearly satisfies 



P 



(N) 



QC 



<P. 



(N) 



(46) 



(47) 



In particular, for pa — Pb = Sp ~ 0, one can show that 
P W _ p (N) 
r QC B 
inequalities [71 Il3| 



Pq V; . Notice that we also have the following 



F- < P (1) < P ( q ] c < F + 



(48) 



where 

P_ := 
and 



l-y/l-F( PA , PB ) 



F + := 



y/F(p A ,p B ) 



F(pa,Pb) := 



Tr 



PAPBy/PA 



(49) 



(50) 



is the fidelity between pa and ps [3. In particular, if 
one of the two states is pure, e.g., p A — Yp) a (vl ? then 
we simply have 



p (i) _ f(Ma (<pUpb) 

^qc- 2 



Formula for Gaussian states 



(51) 



Let us now specialize the problem of quantum state 
discrimination to Gaussian states of n bosonic modes. In 
this case, the quantum Chernoff bound can be expressed 
by a relatively simple formula thanks to the normal mode 
decomposition (x, S, {vk}) of Eq. (Ti"6|) . 

Theorem 2 Let us consider two arbitrary n-mode Gaus- 
sian states pa and ps with normal mode decompositions 
(xA,SA,{afc}) and (x B , S B , {Pk})- Then, we have 

Q s = g s cxp{-id T [V j4 ( S )+V i3 (l- S )]- 1 d} , (52) 

where 



2" EI G s (a fe )G!_ s (/3 fc ) 

Q k=1 
s V det [Va(s)+Vb(1-s)} 

In these formulas d := x^ — x B and 



Va(s) = S A 
Vb(1-«) = Sb 



0A s (a fc )I fc 

:=1 
n 

0Ai- s (A)I 



3 A 1 



,fc=l 



(53) 

(54) 
(55) 



Proof. By applying Lemma Q] to Eq. (|45[) . we get 

Q s =NTr[p A (s)p B (l-s)} , (56) 

where 

n 

N := (Trp s A ) (Trp B - s ) - TJ G s (a k )G^ s ((3 k ) , (57) 

fe=i 

and pa(s), p B (l — s) are two Gaussian states defined ac- 
cording to Eq. ([2"3"| . In particular, the CM's of these 
states are given by Eqs. ([54")) and ([55)1 [according to 
Eq. (|32p ]. For an arbitrary pair of n-mode Gaussian 



5 



states p, p' with characteristic functions x, x' and mo ~ 
ments (V, x) and (V, x'), we have the trace rule 



d 2n i 



R 2 " 



2 n - 



exp 



±(x-x') T (V + V') _1 (x-x') 



^det (V + V) 



(58) 



Then, by using Eq. ([58]) into Eq. (|56|) . we easily get 
Eqs. (J52J) and (JS3J). ■ 

Thanks to the previous theorem, the Chernoff quan- 
tity Q s can be directly computed from the normal mode 
decompositions (x^, Sa, {dk}) and (xg, S B , {/3fc}) of the 
Gaussian states. Notice that this theorem is already 
contained in Ref. 0] but here the formula of Eqs. (j52"|) 
and (|53p is conveniently expressed in terms of the sym- 
plectic spectra {ctk} and {Pk}- 

In applying this theorem, the more difficult task is the 
algebraic computation of the symplectic matrices and 
Sb to be used in Eqs. ([54)) and (|55|) . In fact, while finding 
the symplectic eigenvalues {vk} is relatively easy (since 
they are the degenerate solutions of a 2n-degree polyno- 
mial), finding the diagonalizing symplectic matrix S is 
computationally harder (since it corresponds to the con- 
struction of a symplectic basis [I3|)- For this reason, it is 
very helpful to derive bounds for the minimal error prob- 
ability pW which do not depend on S and, therefore, 
are much easier to compute. 



IV. COMPUTABLE BOUNDS FOR 
DISCRIMINATING GAUSSIAN STATES 

Let us derive bounds that do not depend on the affinc 
symplectic transformations (x^S^) and (xb,Sb), but 
only on the symplectic spectra {afc} and {/3k}. This is 
possible by simplifying the determinant in Eq. (|53p in- 
volving the two positive matrices ~Va(s) and Vg(l — s). 
Such a determinant can be decomposed into a sum of 
determinants by resorting to the Minkowski determinant 
inequality (l8| . In general, such an algebraic theorem is 
valid for non-negative complex matrices in any dimen- 
sion (see, e.g., Appendix [A"]l . In particular, it can be 
specialized to positive real matrices in even dimension 
and, therefore, to correlation matrices. 

Lemma 3 Let us consider a pair of In x In real, sym- 
metric and positive matrices K and L. Then, we have 
the Minkowski determinant inequality 



[det(K + L)] 1/2n > (detK) 



>l/2n 



(det L) 



l/2n 



(59) 



By combining Theorem [5] and Lemma |31 we can prove 
the following 

Theorem 4 Let us consider two arbitrary n-mode Gaus- 
sian states pa and ps with symplectic spectra {cik} and 



{(3k}- Then, we have the "Minkowski bound' 1 

1 N 



QC - 2 



inf M s 

0<s<l 



where 
M s := 4™ 
and 



Y[ * s (a fe ,/3 fe ) + Y[ ¥i-«G9*,a*) 

Lfe=l k=l 



(60) 



■ (61) 



%(x,y) := [*+(aO$i_,,(v)] 



l/n 



(62) 



Proof. By taking the nth power of Eq. ([55]) , we get 



[det(K + L)] 1/2 > (detK) 1/zn + (detL 



l/2n 



>l/2ra 



(63) 



Such inequality can be directly applied to the CM's 
V A (s) and V fl (l - s) of Eqs. (M) and ^B. Then, by 
inserting the result into Eq. (|53[) . we get 



Qs < 



2™ n G s (a fc )G!_ s (/3 fc ) 
k=l 



{[detV A (s)]^ + [det V B (1 



(64) 



By using the binomial expansion and the relations 

n 

detV A (s) = H[A s (a k )} 2 , (65) 

k=l 

n 

det V B (1 - s) = [] [M-s((3k)} 2 , (66) 



k=l 



we get 



' = ^W|A G s (a fe )Gi_ s (/3 fe ) ' (67) 
Now, by using the Eqs. (fT9"j) and ([2H)l . we get 



* /In 



i=0 




n *J(«*)*r-.C9*) 



.k=l 



,k=l 

and, by using Eq. (|62p . we derive 

n 



(68) 




.fc=i 



An Z_/ 

i=0 

JJ*i_,(/3 fc> a fc ) 
fc=i 

n 

JT * s (a fe , /3 fc ) + FJ a fe ) 



.fe=i 



fc=i 



(69) 
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Latter quantity corresponds to the inverse of the one in 
Eq. (|6ip . Now, since every convex combination of posi- 
tive matrices is still positive, we have that exp{- • • } < 1 
in Eq. Then, we get 



< Qs < M s 



(70) 



leading to the final result of Eq. (|60|) . ■ 

The basic algebraic property which has been exploited 
in Theorem 0] is the concavity of the function " 2 \/ det" 
on every convex combination of 2n x 2n positive matrices 
(like the correlation matrices). Such a property is simply 
expressed by the Eq. ([59"| of Lemma [H It enables us to 
decompose the determinant of a sum into a sum of deter- 
minants and, therefore, to derive the bound in the "sum 
form" of Eq. (|6ip . Now, thanks to the Young's inequal- 
ity (l9j . every convex combination of positive numbers 
is lower bounded by a product of their powers, i.e., for 
every k, I > and < 6 < 1, one has 



6k + (1- 0)1 > fcV" 



(71) 



As a consequence, every sum of positive determinants 
can be bounded by a product of determinants. Then, 
by applying the Young's inequality to Theorem [51 we 
can easily derive a weaker bound which is in a "product 
form". This is shown in the following corollary. Notice 
that this bound can be equivalently found by exploiting 
the concavity of the function "log det" on every convex 
combination of positive matrices, (see Appendix [A] for 
details.) 

Corollary 5 Let us consider two arbitrary n-mode 
Gaussian states pa and ps with symplectic spectra {a k } 
and {Pk}- Then, we have the "Young bound" 



M (N) < y(N) ,. 



1 



inf Y s 

0<s<l 



A' 



where 



and 



Y s :=2 n l[T s (a k )T^ s ((3 k ) 



k=l 



T p (x) := [{x + l) 2p - (x - 1) 



2l>] 



(72) 



(73) 



(74) 



Proof. From Eq. (j7"Tj) with 9 = 1/2, we have that every 
pair of real numbers k, I > satisfies 

k + l>2\fkl. (75) 
Then, for positive K and L, we can apply Eq. (|T5[) with 



k := (det K) 1/2,i > , / := (det L) 1/zn > . (76) 



l/2rt 



This leads to the further lower bound 



(detK 



l/2n 



(det L) 



l/2n 



> 2" [detK det L] 



1/4 



(77) 



By applying Eq. ([77]) to the CM's V A (s) and V B (1 - s), 
and inserting the result into Eq. (|64p . we get 



Ms < n 



G s (a k )G^ s (f3 k ) 



k J[ VdctV A (s)dctV B (l-s) 
Then, by using Eqs. (|65|) and (I66p . we can write 

G s (a k )G^ s (l3 k ) 



Y s . (78) 



Y S = U 



(79) 



Exploiting Eqs. (fT9|) and (|20p. we can easily derive 
Eq. d73j) , where 



T p (x) := [*+(x)9-(x)] 



(80) 



also equivalent to Eq. ([74)1 . Finally, since M s < Y s , the 
result of Eq. ([72^1 is straightforward. ■ 

As stated in Theorem 2] and Corollary [5[ the two 
bounds M^ N ' and Y^ N ' depend only on the symplectic 
spectra {a k } and {(3 k } of the input states pa and ps- 
As a consequence, such bounds are useful in discriminat- 
ing Gaussian states which are unitarily inequivalent, i.e., 
such that 



p A f Up B tf 



(81) 



for every unitary U . In fact, since pa and pb are Gaus- 
sian states, every unitary U such that pa — V pbU^ must 
be a Gaussian unitary U = t/x,s- Its action corresponds 
therefore to an affine symplectic trasformation (x, S), 
which cannot change the symplectic spectrum (so that 
{«fc} = {Pk})- Roughly speaking, the previous bounds 
are useful when the main difference between pa and ps 
is due to the noise, whose variations break the equiv- 
alence and are stored in the symplectic spectra. This 
situation is common in several quantum scenarios. For 
instance, in quantum illumination [20l[2l|. where two dif- 
ferent thermal-noise channels must be discriminated, or 
in quantum cloning, when the outputs of an asymmetric 
Gaussian doner [22j must be distinguished. 



A. Discrimination of single mode Gaussian states: 
an example 

Let us compare the bounds of Theorem Q] and Corol- 
lary [5] with the fidelity bounds of Eq. (|48[) in a simple 
case. According to Ref. [23j, the fidelity between two 
single- mode Gaussian states pa and pb, with moments 
(VajX-a) and (Vb,xb), is given by the formula 



F(Pa,Pb) 



where 



exp [-id^V^+VBrM] , 

(82) 



A:=det(VA+V fl ) , 5:= (det V A - l)(det V B - 1) , 

(83) 



7 



and d := -x^. Let us discriminate between the two 
single-mode states: pa = c(l) = |0) (0| (vacuum state) 
and pb = cr(P) (arbitrary thermal state). In such a case, 
it is very easy to compute the infima of M s and Y s in 
Eqs. ([60)) and l|72|) . respectively. In fact, by exploiting 



(84) 

we get 



M, 



(85) 

whose infima are taken at s = and s = 1, respectively. 
As a consequence, for a single copy of the state, we have 

M( 1 ) = (l + / 9)- 1 , F« = ^= . (86) 



At the same time, we have 

F(p A ,p B ) = 2(1 + /3)" 1 , 

which implies 



(87) 



F_ = i-i ^ F 1 

2 2 y /? + 1 ' + ^2(1 + /?) 



By using Eq. (|5"Tj) . we also derive 



(89) 



V. CONCLUSION 

We have considered the general problem of discriminat- 
ing two Gaussian states of n bosonic modes, supposing 
that N copies of the state are provided. To face this prob- 
lem, we have suitably recasted the formula for the quan- 
tum Chernoff bound given in Ref. Q ■ By combining this 
formula with classical algebraic inequalities (Minkowski 
determinant inequality and Young's inequality) we have 
derived easy-to-compute upper bounds whose computa- 
tional hardness is equivalent to finding the symplectic 
eigenvalues of the involved Gaussian states. Since these 
upper bounds depend only on the symplectic spectra, 
they are useful in distinguishing Gaussian states which 
are unitarily inequivalent. This is indeed a common situ- 
ation in various quantum scenarios where the noise is the 
key element to be discriminated. For instance, the dis- 
crimination between two different thermal-noise channels 
is a basic process in quantum sensing and imaging, where 
nearly transparent objects must be detected 13, El- Po- 
tential applications of our results concern also quantum 
cryptography, where the presence of noise is related to 
the presence of a malicious eavesdropper. 
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As evident from Fig. [T] the Young bound is tighter 
than the fidelity bound F + , while the Minkowski bound 

M W coincides exactly with Pqq in this case. 



o 




5 10 15 20 



P 

FIG. 1: Behavior of the various bounds Y w , M (1 \ p££, F+ 
and F- versus the eigenvalue (3 in the discrimination of a 
thermal state o"(/3) from a vacuum state. Notice that = 
Pqc in this example. 



APPENDIX A: BASIC ALGEBRAIC 
INEQUALITIES 

For completeness we report some of the basic algebraic 
tools used in our derivation (see also Refs. [l8l. Izj. [25|). 
Here we denote by M(m, C) the set of m x m matrices 
with complex entries. 

Theorem 6 (Minkowski determinant inequality) 

Let us consider K, L S M(m, C) such that = K > 
and = L > 0. Then 

[det(K + L)]» > (detK)~ + (detL)™ . (Al) 
More generally 

{det [0K+(1 - 0)L]}~ > 9 (det K) ™ +(1 — 9) (det L) ™ , 

(A2) 

for every < 9 < 1 . 

By taking the mth power of Eq. (|A1[) , it is trivial to check 
that 

det(K + L) > detK + dctL . (A3) 

Then, by using the Young's inequality of Eq. (I71| . we 
can easily prove the following 
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Corollary 7 Let us consider K,L G M(m, C) smc/i £/iaf 
Kt = K > and V = L > 0. Then 

det [0K + (1 - 0)L] > (dot K) 9 (det L) 1_e , (A4) 

/or every < < 1 . 

Proof. By setting k := (detK) 1/m > and I := 
(detL) 1 ^" 1 > 0, we can apply Eq. (|7Tj) to the right hand 
side of Eq. (|A"2|) . Then, we get the result of Eq. ljA4| by 
taking the mth power. ■ 

Notice that, by taking the logarithm of Eq. (|A4[) . we get 
/ [0K + (1 - 0)L] > 0/ (K) + (1 - 0)/ (L) , (A5) 



where 



/(M) := logdet(M) . (A6) 



In other words, the Corollary [7] states that the function 
"log det" is concave on convex combinations of positive 
ma trice s. The Theorcm[6]instead states that the function 
" y/ det" is concave on convex combinations of m-square 
non-negative matrices. 
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